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We calculate the fermion Green function and particle-hole susceptibilities for a degenerate two- 
dimensional fermion system with a singular gauge interaction. We show that this is a strong coupling 
problem, with no small parameter other than the fermion spin degeneracy, N. We consider two 
interactions, one arising in the context of the t — J model and the other in the theory of half-filled 
Landau level. For the fermion self energy we show in contrast to previous claims that the qualitative 
behavior found in the leading order of perturbation theory is preserved to all orders in the interaction. 
The susceptibility XQ a^t a general wavevector Q ^ 2pF retains the fermi-liquid form. However the 
2pF susceptibility X2pp either diverges as T — * or remains finite but with nonanalytic wavevector, 
frequency and temperature dependence. We express our results in the language of recently discussed 
scaling theories, give the fixed-point action, and show that at this fixed point the fermion-gauge-field 
interaction is marginal in d = 2, but irrelevant at low energies in d > 2. 



I. INTRODUCTION 

The problem of fermions in two dimensions interact- 
ing with a singular gauge interaction has arisen recently 
in two phiisical contexts. One is the "gauge theoij¥p 
approaches to the t — J model which has been arguecU'cl 
to contain the essential physics of high Tc superconduc- 
tors. Tba, other is the theory of the half- filled Lan- 
dau levelBu. In both cases one is led to the theoreti- 
cal problem of a degenerate Fermi gas interacting with 
a gauge field characterized by the propagator D{uj, k) ^ 
('X' + |fc|^+^)~^. This -notation is conventional; x = 1 
in the t — J model caselll and, because of the unscreened 
Coulomb interaction, x in the 1^—1/2 case consid- 
ered by previous authorsoQ. If the Coulomb interaction 
in two dimensions were screened, e.g. by a metallic gate, 
the model with x = 1 would apply even to the = 1/2 
case. The three dimensional version of this model with 
X = 1 was shown by ReizerQ to describe electrons in met- 
als interacting magnetically via a current-current inter- 
action. The highly singular behavior of the gauge propa- 
gator at small w, k complicates the analysis of the theory 
and has led to conflicting claims in the literature. In this 
paper we present what we believe is a correct treatment 
of the low energy properties of the theory. 

We study fermions interacting with a gauge field a, 
and also with each other via a short range interaction 
W. We assume the simplest form of the intcaraction with 
the gauge field allowed by gauge invarianceOl: 

H = ^ cj;^e(p)Cp„ -I- 4+fc/2,^^kV(p)Cp_fc/2,^ 
p,cr p,k,cr 



+Y1 ^4l,aCp2,a43,'TCp4,<T'5(Xlp0 + Tt/^i' (1) 



where we omitted higher order terms in the gauge field 
a which lead to less infra-red singular effects. Here, as 
usual, /^^ = dfj^tti, — d^a^, go is the bare fermion-gauge 
field interaction constant, a = 1..N is a spin index and 
^ ~ §p' '^^'^ term comes from integrating out high 
energy processes. In the t— J model the spin degeneracy 
N = 2; however, it will be convenient to consider general 
values of N because the limits iV — > and TV — + oo are 
solvable. Indeed, as we shall see N is the only expansion 
paramter of the model. 

In the t — J model the fermion operators c| ^ create 
"spinous" which are chargeless, spin 1/2 quanta. Be- 
cause the spinous have no charge, there is no long-range 
Coulomb term in bF. In this representation the charge 
is carried by different, spinless quanta which obey bose 
statistics ("holons"); we shall not consider their prop- 
erties in this paper. The Hamiltonian (|l]) describes the 
magnetic properties of the "spin-liquid" state oi t — J 
model. For a more detailed discussion of the physical 
situations to which jthis model may apply, see, e.g., the 
review paper by LeeB. 

In the V — \/2 case the spin degeneracy iV = 1 and 
one must add to the Hamiltonian (|^) additional terms 
containing the Coulomb interaction and Chern-Simons 
term. This changes the k^ term in gauge propagator to 
|fc|; the effects of tthis change will be discussed below in 
Section IV. Other authoraJ have found it convenient to 
consider a continously varying exponent |A;|^+^; we find 
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that the behavior for all a; > is the same as that for 
X = 1 except for minor changes in exponents. The x = 
case is exceptional because there a controlled expansion 
for the infrared behavior exists for any N. We give results 
for general x > in Section V. For the rest of this Section 
we explicitly consider only the "spin liquid" case. 

Treating the fermion-gauge field interaction in (|l|) by 
perturbation theory leads immediately to two effects. 
Dressing the gauge propagator by a particle-hole bubble 
leads to the propagator 



D(uj,k) 



1 



Npo\uj\ 



fc2 



(2) 



Here the first term in the denominator is due to Lan- 
dau damping of the gauge field, and po is the curvature 
of the Fermi surface at the point where the normal to 
the fermi surface is perpendicular to k. The second term 
in the denominator has contributions from the f,^„ term 
in the effective action and from the fermion diamagnetic 
susceptibility; in this term we have redefined the inter- 
action constant so that the characteristic energy scale 
remains finite in the limits N oo and N ^ which we 
consider below. 

Using the gauge field propagator to calculate the 
fermion self energy E in the first order of the pertur- 
bation theory one findsllj 



1^0 



1/3 



(3) 



where the energy scale loq is defined in terms of g, po and 
the fermi velocity vp via 
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2\/3 



2vW 



(4) 



For high- Tc materials {N = 2) g^ was estimated to be 
6\/27rm (where m is the fermion mass). This leads to 
ujQ ~ 500 K if the fermion bandwidth is of the order of 
2J. 

The dramatic effects found in the leading order of per- 
turbation theory lead one to question whether the pertur- 
bation theory makes, sense. Several different treatments 
have appearedQ^Efo. The appearance of N in the de- 
nominator of the gauge field propagator (^) suggests that 
the theory should have a tractable N ^ oo limit and 
that a 1/A^ expansion about this limit is well behaved. 
The N ^ oo limit and the leading 1/A^ corrections to 
the ferjpion propagator have bcpa studied by loffe and 
Larkintl, by ReizerQ and by Ledl3 but the higher order 
corrections and the issue of convergence of the expansion 
have not to our knowledge been previously examined. 
We present this analysis in Section II of this paper. We 
find that the 1/N expansion is indeed well defined and 
the leading order results are qualitatively correct for all 
physical quantities except the 2pp susceptibilities, which 
acquire additional non-analytic power law dependences 
with exponents which vanish as ^ oo. 



In order to explain the idea of the analysis we need to 
introduce some notation and establish typical values of 
momenta and energies involved in virtual processes. The 
typical momentum k^^ transferred in a low energy pro- 
cess affecting a fermion with energy lo and momentum p 
is found from the gauge field propagator, eq. (0), to be 



2 \ 1/3 



(5) 



It is convenient to choose Cartesian coordinates in mo- 
mentum space so that fcj^ is the change of the momentum 
of the fermion along the Fermi surface (i.e. perpendicu- 
lar to p) and k\\ perpendicular to the fermi surface (i.e. 
along p). At low energies fcy ~ — Tj{uj)\/vf becomes 
much less than k± which is determined by the gauge field 
propagator (||), i.e. /cj_ ^ k^i- 

Qualitatively the small value of the higher order cor- 
rections at large N can be attributed to a comparatively 
large typical momentum transfer (||) in this limit as fol- 
lows. In a typical virtual process fermion probes only a 
small patch of the Fermi surface of the order of k^^. The 
curvature of this piece of the Fermi surface is important 
if the change in the fermion energy induced in such a 
virtual process {vFk'^/{2po)) is large compared with the 
imaginary part of its self energy (|^). Comparing the two 
we find 



2poS(i)(6) 



N 



(6) 



Thus, in the limit of large N the curvature of the Fermi 
surface becomes important and we expect that the scat- 
tering becomes essentially two diaiensional. In this case 
the usual phase space argumentsEj show that all crossing 
diagrams are small in so that a expansion is 

possible. 

An alternative solvable limit, namely—jV — > 0, was 
pointed out by loffe, Lidsky and Altshulerliil. In this limit 
the curvature of the Fermi surface becomes unimportant 
and the terms proportional to fc^ in the denominators of 
the fermion Green function are negligible. When these 
terms are dropped the Green function does not depend 
on fcj^, which enters only via the propagator of the gauge 
field (g). Thus, in any diagram one can integrate inde- 
pendently all the gauge field propagators over k±. The 
gauge field propagator becomes 



D^^°\uj) = J D{uj,k){dki_) 



(7) 



,1/3 



^i^l^^ is the effective 



where g = {2'Kg^/po 
interaction constant. Note that fcy does not appear be- 
cause it is negligible relative to cu for the reason given 
below eq. (^. 

After these transformations diagrams which do not 
contain fermion loops (except for those loops implicit 
in the gauge-field propagator) become the same as in a 
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ID theory with a retarded interaction given by (|7|) and 
the diagrams which contain loops are negligible. There- 
fore, in the limit — > the theory can be solved by 
bosonization methods. Moreover, by reproducing this 
solution using the diagram technique we find that one 
dimensional results depend crucially on the exact can- 
cellations specific to ID models and that at any N ^ 
these cancellations are not exact. These observations al- 
low us to obtain some information about the behavior at 
<C 1. The analysis of A*" ^ limit is given in Section 
III. The solution at TV ^ turns out to be very simi- 
lar to the results for the fejjmion propagator obtained by 
Khveschenko and StampEil viapoikonal methods and by 
Kwon, Marston and IIoughtoEll3 via a two dimensional 
bosonization. From our results we see that these calcula- 
tions are only valid in the strict ^ limit, so that the 
claim of these authors to have calculated the low energy 
behavior exactly at A^ = 2 or for the half-filled Landau 
level is in disagreement with our results. 

A third theoretical approach involves scaling equations 
constructed by eliminating |high energy degrees of free- 
dom. J. Gan and E. WongEj derived an action for the 
gauge field alone by integrating out the fermion degrees 
of freedom in Hamiltonian and then showed that this 
action has an infra-red stable weak coupling fixed point 
in 2 spatial dimensions. From this they concluded that 
Eq. (g) gives the correct asyipototic form of the gauge 
field propagator. Kwon et alli3 obtained the same re- 
sult via bosonisation. Our results for finite A^ imply that 
"correct asymptotic form" means that the scaling uj ^ 
is preserved, as is the behavior in the limits and 
fc ^ fccj, but not the precise functional form when k ^ k^^. 
An alternative scaling treatment was given by-Ciayak and 
WilczekQ, extending previous work of ShankaJlj on short 
range interactions. Nayak and Wilczek wrote a scaling 
relation for an action based directly on Eq. (|l|). They 
concluded that for the v = \/2 problem in d = 2 the 
fermion gauge field interaction is marginal and in the 
"spin liquid" case it is relevant, so that no statements 
can be made until the strong coupling fixed point is 
found. However, our results imply that the strong cou- 
pling fixed point has a straightforward interpretation: in 
the "spin liquid" case in d spatial dimensions the bare 
scaling e ^ vk ^ vpk'^/{2pQ) is replaced by the new 
scaling e''/^ ~ VFk\\ ~ VFk'j_/{2po) found from the lead- 
ing order gauge field corrections to the fermion propaga- 
tor. In d > 2 we show that any additional corrections 
from the fermion-gauge- field interactions are irrelevant. 
In d = 2 we show that the corrections are marginal at 
X > and lead to new power laws only in the 2pp suscep- 
tibilities. In the case of half-filled Landau level {x = 0) 
these power laws are replaced by a much weaker singu- 
larity. For the case of the half-filled Landau level with 
unscreened Coulomb interaction our results amount to a 
justification of the leading-order approach of Halperin et. 
al.El. The interpretation of our results in terms of scaling 
theory is discussed in Section V. 



Section VI is a conclusion in which the physical inter- 
pretation of our results is discussed. 

After this manuscript was completed we learned of 
two preprints reporting results very similar to some of 
those reported here. Kim, Furusaki, Wen and Ledlll cal- 
culated particle- hole bubbles at small q to order 1/A^^ 
in the spin liquid model, finding, as we did, that the 
fermi liquid forns-js not modified by the gauge interac- 
tion. PolchinksiliS performed a scaling analysis of the 
large-N spin liquid model and concluded, as do we, that 
the curvature of the fermi surface is important and that 
Migdal-type arguments justify the results of the leading 
order perturbation theory calculation. He also obtained 
our result, eq. (|2^), for the renormalization of the 2pF 
component of the four fermion interaction. 

II. LARGE N LIMIT 

This section will show that in the limit N ^ oo the 
leading contribution is given by the diagrams with the 
minimal number of crossings; this will allow us to con- 
struct a perturbative series in 1/N and obtain physical 
results in the leading orders of this expansion. We find 
that to all orders in the expansion the self energy remains 
proportional to e^/'^, that all particle- hole susceptibilities 
except those at |Q| — 2pp retain the usual Fermi liq- 
uid form and that correlators at 2pF momentum transfer 
acquire an anomalous power law dependence. 

In order to develop a consistent large A'^ expansion for 
the Hamiltonian we must take N ^ oo limit so that 
the interaction parameter in (|^) remains constant. At 
N — oo the only diagrams that survive are the RPA 
bubble graphs shown in Fig. |l|. These bubbles screen the 
1/k^ behavior of the gauge field. Because the gauge field 
is transverse, it is not completely screened and the result 
isEq. d). 

We now consider the 1/N corrections to the fermion 
propagator. These are shown diagrammatically in 
Fig. H. The self energy appearing in the leading diagram 
(Fig. ga) was given in Eq. (||). One sees that at energies 
less than cuq or length scales longer than vp/uo, the self 
energy becomes larger than the inverse of the bare Green 
function. We have chosen the way the limit A?^ — > c» is 
taken so that the scale loq remains constant. Because the 
first correction is of the order of 1 and not of 1/A^, care is 
required in carrying the 1/A'^ expansion to higher orders. 

Now consider the 0(1/A^^) terms. The first of these 

(Fig. |b) scales as [S^^H^)] ^ /e ~ e^^^- Direct calculation 
using bare fermion propagators shows that the second 
term (Fig. ||c) scales as e (up to logarithms). Specifi- 
cally: 

^l('^P\\) = 7T^'^''-'''^P\0 fln .^y (8) 
" [2'kYN^ " \ \f. + ivpp\\\j 

where c ~ 3.28 and py = \p\ — pp. This shows that in 
the low energy limit the self energy is more singular than 
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the vertex correction and should be summed first. To 
calculate to higher order in the \/N expansion we should 
therefore use the Green function G^^^ given by 

G«(6,p) = . ^ (9) 

le — vfP — ' (ej 

with the self enegy 'E^^'>{e) is given by (||). In fact this 
G^^^ solves the self-consistent Eliashberg equation 
/ DG also, because E^^^ (e) is momentum independentEj. 
Therefore, the rainbow graphs have been summed and 
we need only to consider graphs with crossed lines such 
as shown in Fig. ||c. 

Returning now to the vertex corrections we reevaluate 
the leading vertex correction, shown in Fig. ||a, using (||) 
for the fermion Green functions. We find that this correc- 
tion is at most of the order of the bare vertex, moreover, 
it is of order (ln[A^]/A^)^ for external momenta of order 
fcij. Explicitly, we find 



-.(li?,) 
- p,0 



VF 



F 



F{x) 
Fix) 



(10) 



0.9422 



a; > 1 
a; = 



Qualitatively, the small value of the vertex correction 
at large N can be attributed to the argument underlying 
the Migdal theorem in the electron-phonon problcmEj, 
namely that the "velocity" of the boson is much less than 
the "velocity" of the electron (by "velocity" in the present 
case we mean However, the argument is more sub- 

tle than in the electron-phonon problem because here we 
have only small angle scattering. To understand how the 
argument goes, consider again the second order crossed 
graph for the self energy Fig. |^, using now ^ for the 
fermion Green function. 



I](2)(e,p) = 4 ^ j G^^\€+UJi,p+ki)G^^\t + UJ2,P+k2)G^^\t + UJi+UJ2,p+ki+k2) 

X D{LJi,ki)D{u;2,k2){(fki(fK2) (11) 
In order to evaluate (O) we integrate over the parallel components of the momenta kin and ^211 obtaining: 



v(2)/ \ _ 2 ' /" DiuJi,kj_i)D{uj2, kj_2)idk±idk±2) 



(12) 



where the prime means that the sum over frequencies is restricted to the region where uji + uj2 + £ has sign opposite 
to + e and uj2 + e and 



(13) 



Clearly, the second order contribution to the self en- 
ergy is at most of the order of 1/N because it contains 
coming from two gauge propagators and N from 
the phase volume (note that k^^ oc N^/"^). In fact, the co- 
efficient of the 1 /N term vanishes because the expression 
under the integral in (^|) is odd in ki± and k2± and the 
leading behavior turns out to be 

(14) 

The reason for the powers of 1/A^ is essentially that the 
phase volume available for the process when all three 
electron lines are on the mass shell is negligible as in 
the usual Migdal argumentsEj, although here the phase 
volume is small only in 1/A^. Note that the non-zero cur- 
vature of the Fermi surface is essential to the argument. 
Note also that in spatial dimension d > 2 the leading self 
energy is e'^^^ so that at any N the small parameter of the 
"Migdal expansion" is e 3 . This is related to the fact, 
to be discussed at greater length in Section V, that the 



interaction is marginal in d = 2 and irrelevant in o? > 2. 
Note that although E^^) j^^g ^j^g 

same form as E(i) in 
the limit Wg^'^jep/'^ 3> vfP\\ it does not have exactly the 
same functional form for Wp'^^jep/"^ ^ vfP\\- 

Thus, at 1 all diagrams can be classified by the 

number of crossings and the sets of diagrams with mini- 
mal number of crossings should be summed^ first, a pro- 
cedure well known from localization theoryllj. The result 
of this summation shows that such diagrams indeed give 
the leading contributions to the higher order terms of 
the perturbation expansion but these contributions are 
not sufficiently singular at low energies and contain ex- 
tra powers of 1/iV. We discuss the calculations leading 
to this conclusion in Appendix B. 

The absence of low energy singularities in the higher 
orders of the perturbation theory implies that the results 
obtained in the leading order are modified only slightly 
by higher order terms. 

The discussion so far has shown that the 1/N expan- 
sion is well defined and has established the qualitative 
form of the fermion propagator. Now we verify that 
higher order corrections in 1/N do not change the qual- 
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itative form of the gauge field propagator. This follows 
from the general considerations of Gan and Wongllj, but 
we believe an explicit derivation is valuable because the 
validity of the approach of Gan and Wong (which in- 
volved integrating out gapless fermions and dealing with 
an action involving the gauge field only) may be ques- 
tioned and because the derivation makes clear that al- 



though the two limits {k^ and k <^ k^) are correctly 
given by Eq. (|^), the precise form for k ^ k^^ is changed 
by higher order diagrams. 

We first consider the leading term D^^^ {oj, q), which 
is obtained by evaluating the polarization bubble (Fig 1) 
but with renormalized fermion propagators. This may 
be written 



1 



• 1/3 1 



;|2/3 



vfPw 



1 



iuj„' \e- 



;|2/3 



VFPW 



(15) 



This may be most easily evaluated by subtracting and 
adding the bare bubble obtained using bare Green func- 
tions in (|l^). In the difference term one may integrate 
over p|| first, then sum over e. The result is 



D^^^{io,q)-D-\LO,q) 



1/3| 

vpq^ 



5/3 



(16) 



Thus, for ^ uJq ^^^{vpq)^^'^, the full propagator is still 
of the bare form (H). Further, we only need this propa- 
gator for UJ - iV-3/2g-2p-lg3 ^ 1/2 ^^^^^3/2^ rpj^^g^ 

the only effect of using renormalized Green functions is 
to reduce the upper frequency cutoff (which enters no 

physically interesting result) from vq to Wg ^^'^{ypq)^^'^- 
A very similar calculation shows that 1 /N vertex correc- 
tion to the polarization bubble shown in Fig.^ is of the 
same order as the self energy correction to the bubble, 
i.e.: 



SD{iu,q) 



1/3| 
PO^O 1^1 

vpq^ 



5/3 



(17) 



Thus, for uj less than the upper cutoff lJq '^^'^ {vpq)'^^'^ the 
renormalization of the gauge field propagator is small. In 
particular, for q k,^ it is smaller than the bare part by 
a factor of order of q. However, the two loop diagram. 
Fig. ||b, leads to a correction which is of the same order 
as the leading diagrams (Fig. |l|) in the infrared limit. We 
do not discuss the details of the evaluation here (except 
to note that the dominant contribution comes when the 
internal gauge-field momentum q is almost parallel to the 
external momentum, fc, and that the fermion loop van- 
ishes if the external frequency w = but are of the order 
of unity if k ^ k^). This diagram therefore does not 
change the scaling or the asymptotic forms in the limits 

3> fcij or fc <C but does change the detailed tj, k de- 
pendence of D{u!,k). This discussion also shows that in 
general the long wave susceptibilities preserve the Fermi 
liquid form for small frequencies. 

So far we have discussed the effects of gauge fields on 
the long wave properties of fermions. Now we turn to the 
effects of gauge field on the fermion vertices with large 
momentum transfer. The corrections to the vertex with 
large but arbitrary momentum transfer |q| pp are gen- 
erally small because of the small phase volume available 



for virtual processes which leave both fermions with mo- 
mentum transfer p -I- q -I- k and p -f k close to the Fermi 
surface. The situation changes only for \q\ close 2pp. In 
this case a virtual process with momentum transfer q 
along the Fermi surface leaves both fermions with mo- 
menta p + q + k and p + k near the Fermi surface. 

The leading contribution in 1 /N to the fermion vertex 
Tq is logarithmically divergent at Q = 2pp; we find that 
higher powers of N contain higher powers of logarithms; 
we sum these logarithms using a renormalization group 
method and find power law singularities in T2pp ■ These 
singularities imply that the calculation of the particle- 
hole susceptibility must be reconsidered. Finally, a sin- 
gular susceptibility near 2pp may be further modified 
by the short range four fermion interaction; therefore we 
must consider also the renormalization of this interaction 
by the gauge fields. 

We begin with the diagrams for Tq shown in Fig. ||. 
The diagrams shown there diverge logarithmically if all 
external momenta are on the Fermi surface, external en- 
ergies are zero and the momentum transfer is exactly 
Q = 2pp. Since the energy only enters the Green func- 
tion via Wg^'^jep/'^, the momentum component across the 
Fermi surface via wj^/cy and the momentum along the 
Fermi surface via k\^ the divergence is cut off by the 
largest among 



(ypk^fl^ 



1/2 



fci/77l3/2 



,1/2 



(18) 



If, say, the largest is the external frequency we evaluate 
the diagrams in Fig. || and get 

^r^,. (c.) = + ^ in3(^)) m (1) r^,^ (19) 

where F^^ is the bare vertex at small scales or large fre- 
quencies. The logarithmic nature of the corrections to 
the effective interaction allows us to sum higher orders 
of the perturbation theory by constructing the renormal- 
ization group equation: 



2pF 



dln(l/w) 



1 

2iV 



1 



2ti^N' 



\^{N) V2pp (20) 
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^Prom 
as 



]) we see that the vertex grows at large scales 



2pp 



-L 2pf 



1 



1 



2N 27r2iV2 



In^(A^) 



O 



1 



(21) 



(22) 



Here we used energy ui for the infra-red cut off assuming 
that it sets the largest scale among lu, e, ey, £±. The re- 
sult (|2^) is derived using a large N expansion. It is also 
of interest to evaluate these diagrams at iV = 2. The 
leading order diagram gives a — 0.25; the sum of the 
diagrams shown in Fig. gb and d gives a = 0.35. 

The power law growth of the vertex at 2pF distin- 
guishes fermions with a gauge interaction from an ordi- 
nary Fermi liquid with short range repulsion and leads to 



anomalous behaviour of the spin correlators at Q = 2pp. 
In the absence of a short range interaction effective at 
2pF (i.e. if the interaction W in eq. |l| vanishes) the spin 
correlator is given by the polarization diagrams shown in 
Fig. ^. The leading contributions in powers of (-^Inw) 
come from the diagrams in which the vertical lines of the 
gauge field do not cross. In these diagrams the leading 
contribution originates from the frequency range (and 
corresponding momentum range, which we have not ex- 
plicitly written) 

eF > t^n > • • • t^l > < • • • < OJ-n < £f 

where lo is external frequency. Therefore, the sum of all 
diagrams is given by the diagram shown in Fig. ^ with 
renormalized vertices (|2l|): 



Ii{LO,q)^ / Gie + u/2,p + q/2)Gie-u/2,p-q/2)[Ti^j^iu;,q)]\dpde) 



(23) 



To evaluate the integral in ( |23|) we note that the main contribution to it comes from the range of momenta and 
energies related by e '--^ e|| ^ e± ^ uj (^8|). Estimating the result by power counting we find that if cr < 1/3 (as occurs 

0, (Z = 2pF. 



for large N) the integral ( |23[ ) converges, but if cr > 1/3 it diverges at ui 
these cases separately and find: 



We evaluate the integral in 



n(w,g) = no 
n(c^,g) = 



Po 



UJQ 



-2a 



+ Cq {\q - 2pf\Io) 



UJQ 



1-3(7 



Cq {\q ~ 2pF\k 



x3ct-1 



cr < 1/3 

C7 > 1/3 



(24) 
(25) 



where the coefficients Cq and are of the order of unity for a curved Fermi surface. Below we shall assume that 
c,j = = 1. Since these coefficients depend strongly on the curvature of the Fermi surface, the case of a fiat Fermi 
surface should be considered separately. We do not discuss it further here. 



The spin polarization bubble (^5|) is equal to the spin 
susceptibility if the effects of the short range interaction 
on the spin correlators at 2pF can be neglected. We jus- 
tify this by showing that a sufficiently weak bare inter- 
action is renormalized to zero by the gauge field. Renor- 
malization of the effective interaction by the gauge field 
occurs in the two competing channels shown in Fig. |^. 
In both channels the corrections diverge logarithmically 
if all external momenta are on the Fermi surface, external 
energies are zero and the momentum transfer is exactly 
Q = 2pF- This divergence is again cut off by the largest 
among e, ey and ej^. If, say, the largest is the external 
frequency we evaluate the diagrams in Fig. ^ and get 



6W{cu)^-2{^-^)\nf^)Wo (26) 



where Wq is the bare interaction at small scales or large 
frequencies. Using the renormalization group to sum 
higher orders we conclude that interaction at 2pF decays 
rapidly at large scales: 



W 



Wo 



(27) 



Certainly, Eq. (^) holds only for sufficiently small bare 
interaction Wq. The decay of the interaction implies that 
the spin susceptibility xi^jl) = n(w,q). For larger bare 
interaction Wq > Wc we expect a transition into an or- 
dered state to occur. We will give the theory of this 
transition in a separate paper. Here we note that at Wc 
the interaction docs not scale and that the basic ingre- 
dients of the theory of the transition are polarizability 
of the fermion system (^,^5|) and the four spin fluctu- 
ation interaction shown in Fig. The renormalization 
of the four spin interaction may be treated in the same 
way as that of H. We find that at large N the leading 
diagrams are those shown in Fig. |[ these lead to a diver- 
gent U with the divergence cut off by the largest among 
|a;/eFp^^, q\\/PF and q±/pF- 



6 



U{uJi,uj2,qi,q2) 



fe) +(^) 



, 6CT + 1 



(28) 



Here we denote — max{\\qi\ — 2pF\\,\\q2\ — 2pf\\), 
Lo = max{\ijJi\, 1^2!) and q±^ = (qi — q2)q where q is the 
unit vector in the direction of qi + q2 . 

To summarize: in the hmit of large N the physi- 
cal properties of the spin liquid resemble conventional 
Fermi liquid with the following important differences: (i) 
the scaling relation between energy and momentum is 
changed to (p^), (ii) spin correlators acquire anomalous 
power behavior ( |24| ) at 2pF, (iii) interaction vertices with 
external field at momentum 2pp are strongly enhanced, 
but (iv) the short ranged interaction between quasipar- 
ticles is suppressed. 



III. SMALL N LIMIT 

In this Section we shall show that in the limit ^ 
the motion of fermions becomes essentially one dimen- 
sional and apply methods borrowed from ID theories to 
obtain physical results which turn out to be qualitatively 
similar to the results obtained in the limit N 00. For 
the fermion propagator the N limit is not singular 
and the ID theory gives qualitatively the same result as 
the 1 2D calculation. We find that for the vertex 

function the N limit is singular because it predicts 
an exponentially divergent vertex function rather than 
the power law derived in the limit N —> 00. We shall 
show that the power law behavior remains valid for all 
finite iV, but the power tends to infinity as — > 0. 

The iV limit is defined via equation (H) for 
D{lu, k). In this formula we take N to zero with 17 con- 
stant. To see why this iV — > limit is essentially one di- 
mensional consider again the second order (in the gauge 
field propagator) contribution to the fermion self energy 
shown in Fig. ^3. Let us perform the integration over 
first. In the limit — > the k± dependence of the 
diagram is controlled by the gauge propagators, which 



implies that the main contribution is at k±^ ~ N^/'^e^/^; 
this means that k\ ~ Ne^/^ is negligible compared to 
the self energy of the fermion (S ~ e^^^). Therefore, at 
A^ — > one may neglect the k± dependence of all Green 
function lines. In addition we may neglect all diagrams 
except those in which all gauge field lines connect two 
electrons moving either nearly parallel or nearly antipar- 
allel. The reason is that if the gauge field couples two 
fermions moving in arbitrary directions all components of 
the transferred momentum are limited by fermion Green 
functions and become small: vp\k\ ~ WQ^'^jep/'^. This 
decreases the phase space volume to fc||A:j^ ^ t^^^ (in- 
stead of e) and decreases the gauge field propagator to 
^ e~'^/^. As a result these processes are small (of relative 
order e^f^) and irrelevant in the infrared limit. 

Very similar arguments apply to diagrams containing 
internal fermion loops other than those contributing to 
the gauge field propagator. Here the reason is that the 
gauge field propagator is small in N^/'^. In the calcula- 
tion of the fermion propagator this smallness was com- 
pensated by the infrared divergence of the integral over 
momentum component perpendicular to the Fermi ve- 
locity. This divergence was cut off by k±^ ~ N'^/'^w^^^ 
cancelling the factor of N^^'^. In the diagrams containg 
fermion loops, momenta in different loops are not exactly 
parallel {k± ~ e^'''^) and integrals over kj_ are cut off at 
kj_ ^ e^/'^; this infra-red cut off does not contain com- 
pensating factors of A^. 

In the remaining diagrams one may integrate the gauge 
field lines over the components of momentum perpen- 
dicular to the Fermi velocity and obtain a one dimen- 
sional theory of electrons (with propagator depending 
on frequency and one component of momentum) coupled 
to the momentum independent but retarded interaction 
D^^{uj) defined in (|). As we shall show below, the re- 
sulting theory can be solved by bosonization methods. 
We shall then use Ward identites to obtain information 
about the behavior at small but non-zero A^. 

Therefore, in the limit A^ — s- the sum of alLdiagrams 
can be found from a mapping to a ID theorytil with the 
action 



S 



VF9\p5,k 



11/3 



(dkduj) 



(29) 



Here pt^x = x'^tx the density operator and a is 
a replica index which runs over K values. We take the 
limit if — > to exclude fermion loops, which we have 
argued to be negligible. In a conventional one dimen- 
sional theory with short range interactions loops would 
be present and would affect the values of the exponents. 
Our A^ ^ limit is defined so that loops are negligible 
in the d = 2 gauge problem. If loops are not negligible 
in the d = 2 gauge problem, their effect is not correctly 
given by the ID theory. 



To compute the fermion propagator we may restrict 
our attention to the right moving particles. The theory 
is then the Tomonaga model with a mtarded interaction, 
and has been solved by bosonizationllil yielding 



G(p,e) = j G{x,t)t 



G{x,t) 



2'k{x — ivpi) I (|x| — isgn{x)vFtY^^ I 



-r(2/3)Zo 



1/3, 
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In the Hmit of low energy and momenta close to the Fermi 
surface G{p, e) acquires a simpler scaling form 



G"°'(E,P) = 



-1 (30, 



vf{p-Pf) \ vI^^{p-pf) 



g{u) = -exp[(-l) 



3/4^3/21 



47r 



Although the Green functions ( pO| ) and (^) have com- 
pletely different analytical structures their qualitative 
properties are similar: both are equal to l/vp\pp — p\ 



in the limit uj, 

,y3u|2/3^, 



1/3, ,2/3 







< v\p ~ pp\ 



l/(cjo^''|ep/3) in the opposite limit lUq' \e\^'^ > v\p 



and both behave as 

l/3|^|2/3 



-PF\ 



We therefore expect a smooth crossover from formula 
( ^ ) to (^) as 0. Both describe overdamped 

fermions with a characteristic energy that scales as (p — 
Pp)^/^. Thus, the limit A^ ^ is not singular fpt the 
fermion Green function. Khveschenko and StampEJ ob- 



tained via eikonal methods a form very similar to (30). 
They claimed their result was asymptotically exact for all 
A^. Our derivation, on the other hand, suggests that the 
precise form depends on two special "one-dimensional" 
features: the neglect of internal loops and the neglect 
of the perpendicular momentum in fermion propagators. 
Both these features are present in the A^ — limit and 
in eikonal approximation of Khveschenko and Stamp, but 
are not present at arbitrary N. Of these two approxima- 
tions the most crucial is the neglect of the perpendicular 
momentum. If the p± dependence of the bare Green 
functions is retained the dressed Green function will not 
have the exponential form (|3^).We do not give the alge- 
bra here but below we apply similar arguments to the 
2pp vertex. If the p± dependence is neglected but loops 
are taken into account, the one dimensional formalism 
will lead to Eq. (^) but with a renormalized argument. 
For this reason we do not believe that the exponential 
form is generic, although the correspondence between the 
A^ ^ and N ^ oo limits leaji.us to believe that the 
scaling e^/^ oc p\\ is. Kwon et a]ll3 also obtained a result 



very similar to ( pO| ) from a two dimensional bosonisation 
method in the problem of half-filled Landau level. Again, 
we do not believe the result is correct for any A > 0. 

We now consider the renormalization of the 2pp vertex. 
In the strictly ID limit A^ ^ all diagrams leading to this 
renormalization coincide with the diagrams of ID Lut- 
tinger model (^9|) which has both right and left movers. 
The Laittinger model can be solved by bosonization. One 
findsU that the renormalized vertex J^p^ grows exponen- 
tially: 



exp 



3.9 



2^|w|V3 



(31) 



In order to understand the reason for such rapid growth 
it is convenient to consider the calculation diagrammat- 
ically. In order to obtain the renormalization of the 2pp 
vertex in a conventional Luttinger liquid with a short 
range interaction vsr one first notices that the first cor- 
rection to the 2pp vertex is logarithmic; then it can 
be proved that renacmalization of vsr cancels with the 
fermion self energycj, so that the leading contribution 
to the 2pp vertex comes from the ladder sum. In each 
block of this ladder one can use bare vertices and Green 
functions; finally, the ladder sum exponentiates leading 
to a power law dependence with exponent determined 
by Vsr. In the present problem the singular interaction 
means that the first correction, ST2pp, to the bare 2pp 
vertex F^^^ is a power law. 



<5F^ 



2pF 



3ff 

27r|a;|i/3 



(32) 



but renormalization of the interaction 



3g 



n , still can- 
eels with the fermion self energy and the series expo- 
nentiates leading to the exponential dependence given in 
( pT| ) . The cancellation of the interaction renormalization 
with the fermion self energy is guaranteed by the Ward 
identity of the ID theory. This relates the exact density 
vertex-Fe^jf (w, g) to the exact Green function G, and 
readso: 



G-\e + iu/2,p + q/2) - G-^{e - uj/2,p- q/2) 



luj — vq 



(33) 



This identity implies that the singular part of the product of the full Green function and the renormalized vertex is 
equal to the singular part of the product of the bare Green function and the bare vertex. 
This cancellation no longer holds in two dimensions. Instead in 2D the Ward identity is 



e,p\ 



G-^{e + Lu/2,p + q/2) - G^^e - w/2,p - g/2) 



(34) 



Here we have distinguished the density vertex F*^ from the two components of the current vertex, and we have written 
the two components of the current vertex in coordinates parallel and perpendicular to vf(p). The gauge field couples 
to fermions via the current vertex. 

In the one dimensional Tomonaga model F"'" is absent and F" = vpV'^ = vpT^^ because the current is proportional 
to the density for fermions moving in one direction. In a general two dimensional theory F°, F" and F"'- are not simply 
related; however, in the present model which has only small angle scattering the identity F" = vpV'^ is still valid up 



to terms of the order of e^/"^ or q\. Further, we show in Appendix A that at sufficiently small q^, F" 
related via: 



and F are 
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r^iuo, q) = B{lo, q)sgn{q^^)T\\ (w, q) (35) 

yy-l/2 2 2 

where a = " ^^^^^ . The range of q over which this result applies is given in Appendix A. 
Using (H) and 1^) in Ward identity (||) we find 

WO). ^ G-'i^ + ^/2,P + g/2) - G-^e - u;/2,p- q/2) 

^ e,p l^' 9j — ■ i — r^T ^ (.J ' j 



The vertex of the two dimensional theory differs from the one dimensional vertex (^3|) by the term proportional to 
q± in the denominator of (^7|). Although this term is small in the limit A'^ ^ 0, it is important because it smears 
the singularity which appears at uj — VFq\\ in the ID theory. ^From ( |3^ ) we can calculate the renormalization of the 
2pF vertex as was done for the strictly ID theory. Consider the diagram shown in Fig. put the external Green 
functions on the mass shell and use (g7|). The result IS 

f 1 1 
ST2pp = / dk\\dk±du!- ; — — — — — -- — — J—. ; (38) 

•' ' II I -Ll V ' ||;i 27T\k±\ + Ari/2g2 F-Ll 



The k integral in ( p8| ) is dominated byfc~fc;^cx|w|"'^/'^; 
for k± in this range we estimate _B(cj,fc|[) ~ which 
implies that the k^\ integral is dominated by tfie region 

<ti k^ while the main contribution to the 

kj_ integral comes from the region k^^ ^ k^. Combining 
these estimates with Eq. (^) gives 



Vp 



5T 



2pF 



1 

N 
1 

N 



diodk I 



(c^fci)l/2 ,^[1 + (fc^/fc^)3] 



duj 



11/3 



3/2 



N 



Inn 



(39) 



For a more precise calculation, including the coefficient 
of see appendix C. These corrections exponenti- 

ate as before leading to a power law form for r2pj!. with 
an exponent a which diverges as A^ ^ 0. Explicitly, we 
find: 



16V2 



97rVA 



0(1) 



(40) 



It is interesting to numerically evaluate the exponent at 
A ^ 2. We find a = 0.56. 

/,From the result for T2pp we may obtain as before an 
expression for the polarization bubble if the short range 
2pF interaction W can be neglected. The calculation of 
W is similar to that leading to Eq. (|27| ) in the previous 
section. One obtains a scaling equation 



6W{iu) = (3{N) In ( ) + 0{W'^ 



(41) 



In a strictly ID theory, (3 = and the leading term in 
the scaling equation is proportional to W^. In our case 
we find for small A^ 



(42) 



PiN) = -\c\VN 



Because the (3 function is negative both at small and at 
large A^ we believe it is negative at any A^. Therefore 
we may again apply the calculation which lead us to Eq. 
( p^ ) for polarization bubble, however, since tr diverges as 
Af ^ , the result is Eq. (|2|). 



IV. HALF-FILLED LANDAU LEVEL 



In this section we treat the singular interaction arguedEl 
to be relevant to the problem of the half filled Landau 
level. The physical problem leads to two new features: a 
Chern-Simons term coming from a singular gauge trans- 
formation which eliminates the explicit dependence on 
magnetic field and a long range Coulomb interaction 
(absent in the spin liquid case begause spinons have no 
charge). In previous treatmentsErQij the Coulomb in- 
teraction was taken to be long ranged. We note that 
in many experimental situations the device containing 
the half filled Landau level may also contain a metallic 
gate which screens Coulomb interaction on length greater 
than a screening length k~^. The resulting gauge field 
propagator which includes the RPA self energy of fermion 
loops and takes into account the dielectric constant e of 
the host semiconductor is 



D{uj,k) 



1 



Pol"l 
2ir|fc| 



k+K 



(43) 



Here u = ^ and the appearance of the 1 / (Svr) instead 
of the conventional 27r may be traced to a 1/ (47r) in the 
coefficient of the Chern-Simons termtl. 

In this section we treat the case k = 0; we expect 
the results to apply if the momenta of interest k'^ are 
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greater than k. In the other hmit, one should use the 
results of the previous section interpolated to iV = 1. 
The momenta are those for which two terms in de- 
nominator of (^3|) are comparable. At temperature T, 
typical frequencies are lu = 2'kT and, if k = 0, we find 
that typical momenta fc^ ~ {STrpokBTi/e^Y^'^ . Using 
a typical Fermi momentum for Ga — Al — As system 
Pq = [inny/^ w 8 X 10^ cm~^ and a typical e = 13 
we find that the unscreened results apply if 

K [cm^i] < 4 X lO^T^/^ (44) 

Thus if at T = 0.1 ii' the screening layer is further than 
1000 A from the 2d electron gas, the unscreened results 
apply. If it is much closer, then one should use the results 
of the previous section interpolated to = 1. 

We turn now to computations using D ( ^ ) with k ~ 0. 



The leading order self energy (Fig. |2|) is 

Here the ellipsis indicates terms which are less singular 
as e ^ 0. Arguments identical to those of section II 
show that E^^) also solves the leading order Eliashberg 
equation, so it sums correctly all rainbow graphs. 

We now argue that higher order crossed diagrams give 
less singular contributions to T,{e,p), so that the lead- 
ing dependence is given exactly by (|4|). Consider the 
leading crossed diagram. Fig. |2p, with the fermion prop- 
agators dressed by the self energy ( ]45| ) . After integration 
over parallel momenta and symmetrization in q_\_i, q±2 
one finds 



with 

A(wi,W2,p|i) = VFP\\ + S^^'(e + tJi + t^2) + S^^'(e + tJi) + £^^^(6 W2) 



The prime on denotes the contraint that sum 

over frequencies is restricted to the region where lo\ -I-W2 + 
e has sign opposite to wi -I- e and L02 -\- e. This constraint 
implies that lo\ and 0^2 cannot vanish simultaneously, so 
no infra-red singularities arise from the frequency inte- 
grals. To extract the infra-red behavior of (^) we may 
replace A by its typical value eln(eF/e) and lo\^i by their 
typical values e. The sum over frequency gives a factor 
of e^. The main contribution to the integrals over fci, 
/c2 is a logarithmic divergence coming from the region 
e < < elne; the final result is 



v(2)/ ^ _ ^'"F ln^ln(ei./e) 



In ej? / e 



(47) 



This is smaller than the leading term by the factor 



In ln(ei?/e) 



ln(ef/e) 



Similar considerations apply to higher order crossed 
graphs. 

Our result, that the leading behaviour at small fre- 
quencies is given exactly by the fitsjt order diagram, is 
reminiscent of the Migdal theoremtj, which states that 
the leading low-frequency behavior of the electron self- 
energy in the electron-phonon problem is given exactly 
by the leading order diagram. The physical fact underly- 
ing Migdal theorem is that the momentum transferred in 
an electron-phonon process is large (of the order oi pp) 
while the energy is small (of the order of Debye frequency 
and much less than vpkp). A very similar argument ap- 
plies here. In the calculations leading to Eq. (E^) the 



energy transferred by the gauge field is small, while the 
integral over momenta is logarithmic and only cut off at 
the scale pp. In the spin liquid case discussed in the 
previous Sections all momentum integrals were confined 
to the region of small momenta. The problem simplified 
only in the large N limit where the range of the momen- 
tum integration became large in A. Thus, the problem 
of half-filled Landau level is analogous to-ihe large N 
limit of the spin liquid case. Kwon et alll3 obtained a 
somewhat different result for the fermion Green function 
via a two dimensional bosonisation method. Their re- 
sult is equivalent to applying our previously discussed 
A^ = bosonisation technique to the half-filled Landau 
level problem. As explained in Section III we do not 
believe this is a correct procedure. 

We now turn to polarization bubble and vertices. As 
in the previously considered spin liquid case, the only 
singularities occur in the 2pp vertices. The leading 2pp 
vertex correction. Fig. |^a, is given after summing over 
parallel momenta by 



(dk) 



1 



Po|e| 
27r|fe| 



87r£ 



|fc| 



|ln< 



Po 



(48) 



Again, the leading contribution to the integral over k± is 
a logarithm coming from the region e < vpk\/pQ < elne. 
Performing this integral and evaluating the sum over fre- 
quencies we get 



2pF 



1 



In" 



In 



■max{T,u},VF{Q - 2ppY /po)^ 



(49) 
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Although it is of only academic interest, we note that 
the higher order corrections may be summed to obtain 
the leading singular behavior. As in the case of the self- 
energy, crossed graphs are less singular than ladder ones. 
As in Section II, the sum of the ladder graphs exponen- 
tiates, leading to 



^2pi, = exp 















V T /. 





(50) 



This weak singularity implies that the polarizability is 
not singular, but the leading frequency and momentum 



dependence is weakly singular. 



V. SCALING 



In this section we recover some of the results obtained 
in previous sections via a scaling analysis. Our principal 
result concerns the properties of the effective action Sef / 
of dressed fermions, coupled to a gauge field a in c? 
spatial dimensions: 



S^ff = I i<h;d'k)¥^J\Luf/^'+-^ - + + I idc.d''q){\q\'+^ + ^)|a(c^, q)\' 

+ g (dwid''fci)(fiw2d''fc2)*J,^_fc^*^2,fc2[vi=^(fci) +vf (fc2)](a(a;i-W2,A:i-/i:2) -H/i.c) (51) 



We find that the fermion-gaugc-ficld interaction g is 
irrelevant for d > 2 and marginal in d = 2. Further, for 
X > in d = 2 the marginality of the interaction leads to 
logarithms only in the 2pp response functions. 

Note that Sef / involves dressed fermions with one- loop 
self-energy E = |ti;|'^/*^2+^^ rather than the linear uj depen- 
dence expected for unrepormalized fermions. As shown 
by Nayak and WilczekB, if the linear lo dependence is 
used in Seff, then the fermion-gauge-field interaction g 
is relevant for d < 2 + x, and is in particular relevant in 
rf = 2 for a; > 0. We argue that the strong-coupling fixed 
point to which the Nayak- Wilczek scaling flows is simply 
the Seff we have written above. The argument has two 
steps. The first is the known resultEj that the first order 
correction to the fermion propagator from the gauge field 
interaction is of the form \uj\'^^^'^^^\ The second step is 
that further corrections do not change the form given by 
the first order correction. We have shown this in previous 
sections by explicit solutions of the model in two limits. 
In this section we give a scaling argument leading to the 
same conclusion. 

We first explain our choice of notation in more detail; 
it comes from the fact, seen in the calculations of the 
previous sections, that a gauge fluctuation of momentum 
q couples primarily to fermions in a patch of the fcrmi 
surface where the fermion velocity vp is perpendicular to 
the direction of q. Therefore, in the effective action we 
have written the momentum dependence of the fermion 
fields using local coordinates defined in a patch centered 
on the point (in d = 2) or strip (in d = 3) of the fermi 
surface where vp is perpendicular to q. In this patch the 
gauge-field-fermion interaction is simplified because the 
transverse component of the gauge field is almost parallel 
to vf , so we may replace the cosine of the angle between 
the gauge field a and vf by unity. 

To see that this construction is reasonable, note that 
from the gauge field propagator in Sef / we learn that at 
frequency uj the important momentum scale is 
^From the fermion propagator we learn that the impor- 



tant momentum scale in the direction perpendicular to 
the fermi surface is \u;\'^/((^+^)- for d > 1 this is always 
much less than the scale defined from the gauge propa- 
gator, so that the momentum transferred from the gauge 
field to the fermion is essentially perpendicular to the 
fermi velocity, and the patch construction is well defined. 
Also from the fermion propagator we see that the k± scale 

d 

is \uj\ , Thus in d > 2 the dependence of the fermion 
propagator on k_i is essential. In d = 2, the momentum 
scale derived from the gauge field and from the fermion 
propagator are the same, and the importance of the cur- 
vature term k^/{2po) is determined by a dimensionless 
parameter (e.g. N). We see that in these arguments 
the curvature of the Fermi surface (specified by Pq^) is 
essential. We shall show below that changes under 
scaling, so one must interpret it as a charge in the renor- 
malization group equations. 

We now discuss the "tree-level" scaling procedure. The 
theory has three coordinates: frequency, fcy , and q (which 
we have shown is the same as k±). All three scale diffec,- 
ently. We choose the scaling of fcy following Shankarllj: 
that is, we imagine integrating out fermions in a shell 
given by A/5 < efc < A about the Fcrmi surface and then 
rescaling the momentum perpendicular to the fcrmi sur- 
face to restore the upper cut off. We then choose the 
scaling of frequency to keep the |a;|''/(2+^) term in the 
fermion action invariant, and then choose the scaling of 
q (which is the same as that of k±) to keep the gauge 
field propagator invariant. Finally, we choose the scaling 
of the fields to compensate for the scaling of the coordi- 
nates and integrals, so that the quadratic terms remain 
invariant. Note that we must interpret J d'^qva the gauge 
field term or the fermion field term as / dgyd'^^^gj^. This 
implies 

3d+l+x 

\If ^ vfrfo 2d 

a ab 1^ 
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Combining all factors we get the following tree-level scal- 
ing equations for the charges 1/po a-nd g: 

d\nb^[ d 



dg 

dlnb 







Therefore, the gauge-field fermion coupling is marginal 
for d> 2 but in c? > 2 the effective curvature of the fermi 
surface grows. For large curva±ure the usual arguments 
leading to the Migdal theoremllj imply that the crossed 
graphs may be neglected, so we may restrict ourselves to 
the leading order of perturbation theory, which gives the 
self energy jcjI^T^. 

Alternatively, one may consider a scaling procedure 
which preserves the form of the fermion propagator. In 
this case one must scale kj_ as 5"^^^ and in c? > 2 both the 
coefficient of the g^"'"^ term in the gauge propagator and 
gauge-field-fermion coupling g scale to zero (indeed the 
tree-level scaling equation for g becomes -^^j = so 
that a manifestly weak coupling fixed point is obtained). 

In d = 2, however, all charges are marginal for all a; > 
and further analysis beyond tree level is needed to de- 
termine which physical quantities are renormalized. For 
our purposes the most efficient method of deriving the 
one-loop renormalization group equations is to use the 
technique of differentiating the one-loop diagrams with 
respect to the upper cut off. The calculations presented 
in the previous sections can be carried over directly to 
show that the only quantities which are renormalized are 
the 2pF vertex T2pp and polarizability Il{2pp). In partic- 
ular, neither po nor g scales in d = 2. Rewriting Eq. ( |l9| ) 
in the notations of this section (here we normalize to fc|| 
and in the previous section we normalized to frequency) 
we find 



dT 



2pf 



d\nb 



-ctF 



2pf 



(52) 



where cr is a number which depends on the fixed point 
values of po and g. Our results of the previous sections 
may be viewed as calculations of the fixed pont values of 
Po and g in the large and small N limits. 

Although there are no logarithmic corrections to the 
fermion propagator, the finite renormalizations gener- 
ated by marginally irrelevant operators do mean that the 
fermion propagator G'~^(e,p|| ,pj^) is not precisely given 

by the form — vp{p\\ +p-^/(2po)) written in Eq. ( |5l| ) 

when |e|'+^, vpp\\ and p^/(2po) are of the same order, 
although the limits when one argument is much larger 
than the others are correctly given. 

Finally, we consider d = 2, x = Q. Here at tree level we 
would conclude that for the action with inverse fermion 



Vpkn 



2po J- 



fc^) and for inverse gauge 



propagator (w 

propagator (|gl , |_^| 
marginal. However, caution should be excercised in de- 
riving renormalization group equations beyond tree level 



— ) the fermion-gauge-field coupling is 



because in the two-loop calculations presented in the pre- 
vious section no terms of order (In A)^ were found so that 
the logarithms found in one-loop order do not sum to 
powers. Instead, the calculations presented in section IV 
show that the asymptotic form of the fermion propaga- 
tor is ja;ln|(jj| — — ^^1: and the 2pF vertices are 
extremely weakly singular (~ exp ^ In^ In 



VI. CONCLUSION 

We have presented a discussion of the low energy prop- 
erties of a system of fermions in spatial dimension d cou- 
pled via a singular gauge interaction with propagator 
D{Lj,q) = {\uj\/\q\ + \q\^+^y\ We found that the fermion 

lifetime scales as jwl^T^ (in 3 > d > 2, a; > 0) and that 
in d = 2 the 2pF fermion polarizability n(a;, q) was non- 
analytic and possibly divergent a,s Q —f 2pF and lo 0. 
Whether or not the susceptibility is divergent depends on 
the value of an exponent, a, which we could calculate only 
in certain unphysical limits. In the spin liquid case x — 1 
extrapolation of our calculated a to the physical limit of 
spin degeneracy N = 2 from two sides yielded estimates 
for a bracketing the critical value CTc = 1/3 above which 
n diverges. In the v = 1/2 case one must distinguish 
between screened and unscreened Coulomb interactions. 
In the unscreened = 0, the self energy is luIhuj 

while the nonanalyticity in the 2pF vertex is very weak: 
exp (i ln^[ln[w]]) and the polarizibility does not diverge. 
In the screened case the results for x — 1 apply with spin 
degeneracy N = 1 and our estimates suggest that the 
2pF polarizability diverges. 

There is a simple physical interpretation for the non- 
analyticities at 2pF' a moving fermion emits a gauge 
field which relaxes so slowly that if at a later time the 
fermion is scattered backwards it meets the gauge field 
again and is able to lower its energy. It is remarkable 
that this physics can lead to an actual divergence of the 
2pF susceptibility if the fermion-gauge-field interaction 
is strong enough. The form of the divergence is given in 
Eq. (|2^ , |25| ) and is controlled by an exponent a which 
can a priori take any value. However, we note that if 

> 7/6, then ^H(u;,q) is infrared divergent. Such a 
divergence is not possible; for example in a magnetic sys- 
tem this would imply that the expectation value of the 
square of the local spin density (Sf) diverges. Therefore 
we believe that for cr > 7/6 some other physics beyond 
the scope of our calculations must intervene to cut off 
the divergence. One mechanism for this feedback can 
be seen in the spin-fluctuation contribution to the elec- 
tron self-energy. For cr > 7/6 this diverges, implying a 
smearing of the Fermi surface which would suppress the 
Fermi surface singularities we have found. However, for 
7/6 > a > 1/3 we believe this critical phase is stable. 

In order to understand the physical properties of the 
critical phase, consider first a translation-invariant elec- 
tron gas (as is realized in the half-filled Landau level). 
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Then Eq. ( p5| ) would predict that the susceptibility di- 
verges as T ^ on a ring of radius 2pF- For fermions 
on a lattice, the situation is more complicated for rea=. 
sons very similar to those analysed by Littlewood et alEll 
in a study of 2pp singularities in a marginal Fermi liq- 
uid picture. First, intead of a circle of radius 2pp one 
obtains one or more curves traced out by the vectors Q 
connecting points with parallel tangents. Second, the 
amplitude (but not the exponent) of the divergent term 
in X2pF varies around the curve due to the variation of 
vp and Pq around the Fermi surface. Third, one obtains 
additional families of curves on which x diverges. These 
are generated by Q -I- G where G is any vector of the re- 
ciprocal lattice. As a result one gets additional peaking 
when members of different families intersect. The result, 
for band structures appropriate to high-Tc superconduc- 
tors will be a susceptibility strongly peaked at particular 
points in g-space which might be qualitatively consistent 
with neutron data for La2^xS'''xCuO,^. In addition, the 
divergent spin fluctuations imply that the Cu NMR rate 
1/TiT oc 1/T2'^-i/3, so the 1/TiT diverges as T ^ 
even at the borderline value a — 1/3. The value cr = 2/3 
would lead to 1/TiT cx 1/T consistent with Cu NMR ex- 
periment on high-Tc superconductors. Of course, if these 



wavevectors where x"{'-^tQ) is maximal are too far from 
the commensurate wavevector (tt, tt), the oxygen l^TiT 
will also diverge, in disagreement with experimentE2l. 

In the half-filled Landau level case with screened 
Coulomb interaction the divergence in the 2pp suscepti- 
bility could in principle be observed in sound propagation 
experiments in which the phonon wavevector is tuned to 
2pp. The divergence should lead to a large damping of 
the phonon which increases as T is decreased. The ef- 
fect should be observable for temperatures and phonon 
frequencies less than a scale luq which we calculate from 
Eqs. (§,§,|3|). We rewrite the expression for ojq in terms 
of the Fermi energy ep = h'^p^/(2m), Coulomb parame- 

2 1/2 T 

ter Ec — ^-^ — and the screening length obtaining 
hiuo « 0.15^ (53) 

Assuming typical numbers for GaALAs inversion lay- 
ers m = O.OTtog, e = 13 and n = 10^^ cm~^ we have 
Ep 50 K and E'^ ~ 40 K so hujo [K] w 10^. Thus 
if the screening length is not too much greater than the 
interparticle spacing, the effect should be observable. 



APPENDIX A: VERTEX AT LOW MOMENTUM 
TRANSFER 



Here we use the Ward identity to derive the exact form 
of the renormalized fermion-gauge-field vertex at low mo- 
mentum transfer q± and g|| (the exact conditions under 
which this form applies will be obtained below). In this 

iujTlpiuj, q) - (?|iF|p(tj, q) - (J_lF^p(w, q) = G 



limit the renormalized vertex becomes singular, and our 
goal is to find the form of this singularity. The expres- 
sion that we shall find is correct at any N for sufficiently 
small transferred momenta q±, q^^ . The value of N deter- 
mines only the range of q±, q\\ over which the expression 
obtained in the limit of very low momenta remains valid. 

The fundamental Ward identity was given in Eq. (jsj); 
we repeat it here for convenience. 

1 (e + c^/2, p + q/2) - Q-^e ^ w/2, p - q/2) 



Here F° is the density vertex and F" and F^ are the 
components of the current vertex parallel and perpendic- 
ular to vf(p). We wish to obtain from this an equation 
relating F° to the fermion Green function. As noted pre- 
viously, in the present model at small w the current ver- 
tex F'I is related to the density vertex F" by F" = vpT'^ . 
We now derive the relation between F-'- and F°. Con- 
sider a high order diagram for F-*" of the type shown in 
Fig. in which one of the gauge field lines connecting 
one external fermion leg to the other is isolated, i.e. not 
crossed by any other gauge field line connecting one ex- 
ternal fermion leg to another. The analytical expression 
has the general form 

- - / Y{id^k,)Y,k^,\{D{k,)Y{G (Al) 

where index i runs over n values corresponding to the 
gauge field lines connecting different legs and we did not 



explicitly write the arguments of the fermion Green func- 
tions G. Label the momentum of an "uncrossed" gauge 
field line by k^, and pick out the term in the sum pro- 
portional to k^a- We show below that all other terms in 
the sum are small. 

In the limit q_\_ Q the fermion Green functions in 
diagrams such as Fig. ^ depend only on the combina- 
tion (fc||a + k\^/2pQ), so their dependence on k_\_a can 
be completely eliminated by the shift fcy fcy — k\/2pQ. 
(Recall that the fcy dependence of D is negligible always). 
After this transformation the only remaining dependence 
on fc_La is in D{k±a)- The remaining integration over k± 
is straightforward: 

J D{k^a)k±adk^a = (A2) 

This integral converges poorly at large k±, because the 
integrand decreases as ~ l/k± at large k±, but is zero 
at q± = 0, because the integrand is an odd function of 
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k±. At any qj_ y^ the dependence of the Green fuc- 
tions on the momentum k± can no longer be neglected. 
Since Green functions depend only on the product q±k±, 
their k± dependence becomes significant only at large 
k± ~ S(e)/(7j^. At smaller k± the dependence of the 
Green function on kj_ can be neglected and the contribu- 
tions from positive and negative k_i cancel each other. 
Because the main contribution to this diagram comes 
from large k±a ^ the Migdal arguments of Sec- 

tions II and III show that at large k± all diagrams in 
which other lines cross the line with large momentum 
transfer are small. 

Now consider any arbitrary diagram for . The corre- 
sponding analytical expression will be of the form shown 
in eq. (|Al| ). Pick out one term in the ^J-i. The di- 



agram will be important only if the gauge field line carry- 
ing this momentum is " uncrossed" in the sense discussed 
above. This justifies the assumption made above that the 
term in the sum proportional to k±a corresponds to an 
"uncrossed" line. Therefore, the diagrams which give the 
dominant contribution at small q± can be represented as 
the diagram shown in Fig. |^b. Here the two blocks in- 
volve gauge field lines which cross each other and the 
double wavy line represents D{k±a)k±a- Since the in- 
tegral over k± in the double wavy line is dominated by 
large k±, the frequency dependence of it can be neglected 
while the dependence on fcy can be neglected always. For 
this reason, the outer block is simply with the bare ver- 
tex F''. The inner block can be also expressed in terms 
of the vertex F°. After some manipulation we find: 



Po 

c-H^, p+k 



(A3) 



{n,q)G{e+uj + n/2,p+k + q/2)G{e+uj-n/2,p+k-q/2) 



[(Pkduj) 



In this equation the vector character of the vertex is expressed via the factor of /cj^ which may be positive or negative. 
Together with Ward identity Eq. (jsj), Eq. (A3) forms a closed system of equations for the vertex r°. To solve it we 
introduce the notation 



Ti-jLU,q) = VFB{q»,q^)Tl^{LU,q) 



(A4) 



Here we have suppressed the dependence of B on the frequency lo, because u; is a dummy variable in the analysis that 
follows. We use the Ward identity to express F'^ in Eq. (^3) through B{qii,qj_), finding: 



g^N^/^ f {(fkdw) G{e+uj + n/2,p+k + q/2) - G(e+cj-0/2,p+fc-g/2) 



PO J m-vp{q^\ + ^)-qj^VFB{q+^,q^) 

Integrating over fcii and uj and scaling the k± variable (fcj^ k±po/q±) we find: 



(A5) 



B{q\\,q±) = a sgnqj_ 



P 



dk 



in 



2Trk in - (g|| + fc) - B{q\\ + fc, q^)q^ 



(A6) 



where f2 = n/vp and 



a 



vpg 



2^1/2 



Po 



(A7) 



is a dimensionless parameter of the order of N^^^. It is 
convenient to consider positive and negative f2 separately. 
The considerations are simil ar s o we consider explicitly 
only r2 > here. Equation (A6) can be simplified if we 
assume that the denominator in it has poles only in the 
upper half plane. We shall show this assumption is self- 
consistent. In this case the integration contour can be 
closed in the lower half plane and the integral equation 
simplifies to the algebraic equation 



1 



2 in - (gjj) - q_LB{qii,qj_) 



(A8) 



Solving it we find 



B{q\\,qi 



in - vpqw + \J{in - vpqwY - 2a\qi_\\n\ 



2qj 



(A9) 



Restoring the notations of Section III and using wf9|| 3> 
n we get the equation (36) announced in Section III. 

Combining the Ward identity ( ^^ with Eqs. (A4) and 
(^ ) we get the final expression for the vertex at low mo- 
mentum transfer: 



G{e + n/2,p+q/2) - G{e-n/2,p-q/2) 



VFq\\ + ^/{vFq\\] 



2Q!|f2||(7j_| — invpq^^ 
(AlO) 



Thus, F^ is substantially enhanced at low momenta 
vpqii « |c^o^'||w|2/3, qj_ < {vFqii)V{a\n\). This range 
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of momenta becomes small at iV ^ 1 because a ^ N^^"^ 
and does not contribute much to the self-energy. 

Th e essential ingredient in the derivation of (|3^) and 
( |A1C| ) was assumption that the momentum fcj^ is suffi- 
ciently large so that crossing diagrams can be neglected. 
We also assumed that k± :§> which allowed us to ne- 
glect the frequency dependence of the gauge field propa- 
gator. Both these conditions are satisfied if q± ^ j^k^. 
In the limit ^ 1 this co nditio n limits drastically the 
range of momentum where ( A10| ) can be applied. 



APPENDIX B: HIGHER ORDER DIAGRAMS IN 

1/iV 

The calculations of Appendix A show that the gauge- 
field-fermion vertex is enhanced at very low momen- 
tum transfer as is evident from Eq. (AlO), This equa- 
tion, however, was derived under the assumption that 
q± ^ jfk^. At larger momentum transfer q± t he cor- 
rections to the bare vertex are small, Eq. ( |A10| ) is not 
valid, instead, the leading corrections are given by the 
first crossing diagram shown in Fig. ^. The straight- 
forward calculation gives Eq. (|lO| ) . The equation ( p^ ) 
crosses over to the equation ( |A10D at q± ^ j^ki^. So, at 



large N the momentum range where the whole series of 
diagrams should be summed is small in morreover, 
this momentum range turns out to be so small that it 
does not contribute even to sublcading order in 1/iV for 
most quantities. For instance, the contribution of this 
range to the self energy is of the order of -j^, whereas 
the leading term which comes from larger momenta is 

of the order of (Eq. (|l|)). Thus, in order to 

obtain the subleading terms of the order of it is 

sufficient to keep only the first crossing diagrams in pho- 
ton propagator. However, in order to obta in all terms of 
the order of one needs to use the Eq. ( AlO ) and the 
crossover formulas (which we did not derive) in the range 

q± ^ jjkij. 

Similarly, in the calculation of the exponent of the 2pF 
vertex, the enhancement of the gauge-field-fermion ver- 
tex at q^"^ j^k^ leads to corrections of the order of 
to the exponent. As we shall see below, the leading terms 
are larger by factors of ln[A^], so in the subleading term 
we again can keep only the simplest crossing diagrams 
shown in Fig. ^ Only the diagrams shown in Fig. ^ 
and Fig. |5|d have contributions which contain powers of 
ln[iV]. Consider the diagram shown in Fig. Its ana- 
lytical expression reads 



5iT2pp (w) = v% / D{r], k)D{n, q)G{T], k)G{n + v,q + k)G{n, q)G{n + uj, -q){d^qd^kdndr]) 



(Bl) 



The contribution containing logarithms of N comes from the momentum range fc,, > > \q±\ > -^k^. In this 

range the self energy parts of the Green functions can be neglected. We integrate over parallel components of momenta 
kn and qn and symmetrize the resulting expression obtaining: 



S1T2PJ, (w) 



dfl I dr] 



dk±dq± 
"'0 ^l?! ^ Q± 



Evaluating the remaining integrals with logarithmic accuracy we get 



1 \n^[N] 



In I - 

LU 



87r2 iV2 

Evaluation of the diagram shown in Fig. ^ is very similar It has analytical expression 

62T2pA^)^4 J D{r^,k)Din,q) 

X G{n, k)G{n + V,q + k)G{n + ri + uj,-{q + k))G{n + UJ, -q){d^qd^kdndr]) 
We integrate over parallel components of momenta, obtaining after symmetrization: 



(B2) 



(B3) 



(B4) 



dfl / dr] 



/3|r>|2/3|ry -t- »p/3 - ql{kj_ + q^fsgn[n{n + 77)] 
(/3|17|4/3 + qi){p\n + r,\^ + {q^ + fc^)4) 

X D{Vl,q±_)D{'q,k±^)dk^dq±_ 



(B5) 



Here (3 = A (^f^ j ■ The first term in the numera- 
tor of this integral is logarithmically divergent; the main 
contribution to the integral comes from the frequency 



range 77 > f2 > w and results in a 2^1ii^(l/w) contri- 
bution which one expects from general renormalization 
group arguments. The second term in the numerator 
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has no contribution from this frequency range due to 
sgn[n{n + 77)], so it results in only one power of ln(l/w), 
instead it contains ln[A'^] coming from the momentum 



range fc^ > |fc_L| > |g_L| > 



In this momentum 



range we neglect the self energy parts of the Green func- 
tions and perform integrals with logarithmic accuracy ob- 
taining: 



S2^2pF (w) 



1 \n^[N] 
47r2 iV2 



In 



(B6) 



Adding the contributions from the diagrams in Fig. Ijb 
(which come with a factor of two) and Fig. ||d we get Eq. 



(0 



APPENDIX C: EXPONENT OF THE 2Pf 
VERTEX IN THE SMALL N LIMIT 

In order to find the exponent of the 2pp vertex in 
the small TV limit we evaluate the first correction to the 
2pF vertex shown in Fig. ||a using the exact gauge-field- 
fermion vertices and then exponentiate the result. This 
prescription is known to work in ID Lutinger model and 
it gives the leading terms of the 1/N expansion. The 
analytical expression for Fig. ||a is 



ST2pj,{uj) ^ ~Vp I G(r/,g)G(ry + w,-q)i:>(?7,g)r,,,,(-,-)r_^,_g(-,-)(d qdr]) 



(CI) 



This expression simplifies if the external f ree 
because in this case G~^{e,p) = in Eq 



uency of the fermion is zero and its momentum is on the Fermi surface 
for the vertex: 



-I- y/{vFq\\Y - 2ar]qA 



rD{r],qj_) 



(C2) 



Here we replaced the exact dependence on the external 
frequency uj by an approximate cut off which is sufficient 
for logarithmic accuracy. Evaluating this integral we get 



where a is given by (p 



a\n I - 



(C3) 



1 L. B. loffe and A. I. Larkin, Phys. Rev. B39, 8988 (1989). 
^ G. Baskaran, Z. Zou and P. W. Anderson, Sol. St. Comm. 

63, 973 (1987). 
^ P. W. Anderson, Science, 256, 1526 (1992). 

* F. C. Zhang and T. M. Rice, Phys. Rev. B37, 3759 (1988). 
^ B. Halperin, P. A. Lee and N. Read, Phys. Rev. B 47, 7312 

(1993) 

® V. Kalmeyer and S.-C. Zhang, Phys. Rev. B46, 9889 (1992). 
M. Reizer, Phys. Rev. 40, 11571 (1989). 

* P. A. Lee, p. 96 in High Temperature Superconduc- 
tivity: Proceedings, eds. K. S. Bedell, D. Coffey, D. E. 
Meltzer, D. Pines, and J. R. Schreiffer, Addison- Wesley 
(Reading, CA: 1990). 

^ C. Nayak, F. Wilczek, unpublished. 
1" P. A. Lee, Phys. Rev. Lett. 63, 680 (1989). 
" L. B. loffe, D. Lidsky and B. L. Altshuler (unpublished). 

H. J. Kwon, A. Houghton and J. B. Marston (unpublished). 
" J. Can and E. Wong, Phys. Rev. Lett. 71, 4226 (1993). 



' D. V. Khveshchenko and P. C. E. Stamp, Phys. Rev. Lett. 
71, 2118 (1993). 

' A. B. Migdal, Sov. Phys. J.E.T.P. 7 333 (1957) and A. A. 
Abrikosov, L. P. Gorkov and I. E. Dzyaloshinsky, Meth- 
ods of Quantum Field Theory in Statistical Physics, 

Dover (1975), section 21, Chapter 4. 

' R. Shankar, Physica A 177, 530 (1991) and Rev. Mod. 
Phys. 66, 129 (1994). 

Y. B. Kim, A. Fumsaki, X. G. Wen and P. A. Lee, unpub- 
lished. 

' J. Polchinski, unpublished. 
' B. L. Altshuler, A. G. Aronov, 
interactions in disordered 

and M. PoUak, North-Holland (1985). 

' I. E. Dzyaloshinsky and A. I. Larkin, Sov. Phys. JETP, 38, 
202 (1974). 

P. B. Littlewood, J. Zaanen, G. Aeppli and H. Monien, 
Phys. Rev. B 48, 487 (1993). 

' A. J. Millis, p. 198 in High Temperature Supercon- 
ductivity: Proceedings, eds. K. S. Bedell, D. Coffey, D. 
E. Meltzer, D. Pines, and J. R. Schreiffer, Addison- Wesley 
(Reading, CA: 1990). 



p. 1 in Electron-Electron 

systems, eds. A. L. Efros 



16 



+ ... = VWVA/ 



Figl 



FIG. 1. RPA sum of bubble diagrams leading to dressed gauge field propagator (denoted by thick wavy line). The solid 
lines with arrows denote fermion propagators and the heavy dots denote the bare gauge field propagator. Whether the fermion 
propagators are bare or renormalized does not affect the result of the calculation. 



a b 



Pig 2 



FIG. 2. Fermion self energy diagrams. The wavy line denotes the gauge field propagator (2) and the solid line the fermion 
propagator. 
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Fig3 



FIG. 3. Correction to fermion-gauge field vertex. The wavy line denotes the gauge field propagator (2) and the solid line the 
fermion propagator. 



<x> 



Fig 4 



FIG. 4. Correction to fermion polarizability. The wavy line denotes the gauge field propagator (2) and the solid line the 
fermion propagator. 
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■< = <] + 



Figs 



FIG. 5. Ladder sum giving renormalization of fermion 2pF vertex T2pf (shaded triangle), a. Leading order in ln[N]/N b. 
and c. Subleading order in l/N. d. The leading order in {ln[N]/N)^ . The heavy dot indicates the bare 2pF vertex, the wavy 
line denotes the gauge field propagator (2) and the solid line the fermion propagator. 



= + <I> + <!l> + 



Fig 6 



FIG. 6. Ladder sums giving renormalization of fermion 2pF polarizability. Notation is the same as in Fig. 5 
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Fig 7 



FIG. 7. Diagrams leading to renormalization group equation for short range fermion vertex W (heavy square). The wavy 
Une denotes the gauge field propagator (2) and the solid line the fermion propagator. 



Fig8 



FIG. 8. Diagram leading to a singular interaction between spin waves. Solid lines denote fermion propagators, shaded 
triangles T2pp vertex. 
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Fig9 



FIG. 9. a. Typical high order diagram for F''- . b. Sum of these diagrams. Broad wavy line represents propagator k±D{LJ, k±) 
with large momentum transfer in perpendicular direction. 



21 



